The variable phase method used to calculate and correct scattering lengths 
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It is shown that the scattering length can be obtained by solving a Riccati equation derived from 
variable phase theory. Two methods of solving it are presented. The equation is used to predict 
how long range interactions influence the scattering length, and upper and lower bounds on the 
scattering length are determined. The predictions are compared with others and it is shown how 
they may be obtained from secular perturbation theory. 

PACS numbers: 



I. INTRODUCTION 

It is well established that knowledge of the scatter- 
ing length of a pair of colliding atoms is important in 
the interpretation of observations of the behaviour of 
cold trapped ensembles of such atoms Q, 0, • Scatter- 
ing lengths are usually calculated by solving the single- 
channel radial Schrodinger equation to find the phase 
shifts for several small values of the wavenumber of rel- 
ative motion and extrapolating from effective range ex- 
pansions 0, or by solving the zero-energy Schrodinger 
equation and evaluating quadratures |5|. 

The Schrodinger equation is a second order differen- 
tial equation. We shall show another method of eval- 
uating the scattering length in which we solve the first 
order equation, a Riccati equation, of the variable phase 
approach to potential scattering [6j. Solving the equa- 
tion is not simple and we shall suggest and illustrate two 
methods that circumvent the inherent difficulties. A dif- 
ferential equation can be obtained for the effective range 
@ but it is not amenable to numerical solution. 

In obtaining the scattering length from the solution of 
the Schrodinger or Riccati equation, we must obtain so- 
lutions at infinite values of the interatomic separation, R. 
In practice we stop the solution at some value i? c , but 
i? c must be chosen very large, and such choices increase 
computation time and accumulated error. We shall show 
how the Riccati equation can be used to predict correc- 
tions, to be applied to the calculated scattering length, 
that compensate for stopping the calculation at finite dis- 
tance i? c ; such corrections enable us to make calculations 
with smaller values of R c while maintaining desired ac- 
curacy 0,0I3- 

The solution of a low energy Schrodinger equation may 
be corrected for the influence of a long range interac- 
tion at separations exceeding R c by perturbation theory. 
However the perturbation expansion is made in terms of 
the ratio of the potential energy to the kinetic energy 
and becomes invalid as the energy of the collision is re- 
duced (as is necessary in the usual methods to obtain a 
scattering length) . This problem may be avoided by use 
of secular perturbation theory . We shall discuss our 
predicted corrections in the light of secular perturbation 



theory and also compare them with correction formulae 
that have been obtained elsewhere [Tol ITTl IT^ . 



II. VARIABLE PHASE THEORY 

The s-wave Schrodinger equation for a pair of colliding 
atom is: 



AR Z 



V(R) 



y(k;R)=0, 



(1) 



where V{R) = 2^V(R)/h 2 , V(R) being the interac- 
tion potential, [i the reduced mass, fi the rationalised 
Planck's constant, k the asymptotic (for large separa- 
tion, R) wavenumber of relative motion and R~ 1 y(k; R) 
the wavefunction. The phase shift 5k is obtained from 
the wavefunction at large R: 



y{k;R) ~ sin (kR + 5 k ) 



(2) 



Suppose that Sk(R c ) is the phase shift appropriate to the 
truncated potential V(R)H(R C — R) where H denotes 
Heaviside's unit step function. The phase shift Sk is the 
limit of 8k(R c ) as R c — > oo. The phase shift Sk(R c ) is 
the phase shift that is determined numerically when the 
potential is assumed to be negligible at R > R c ; it has 
long been studied, is called the variable phase function 
and satisfies a first order differential equation in R c 6] . 

The potential for a pair of scattering alkali atoms is 
dominated at long range by the van der Waals interaction 
—C§R~ 6 and the well known problem of finding ultra- 
low energy phase shifts from Eqs. (JIJ and (0) is that R c 
must be chosen very large to ensure that —CqR~ Ci has 
negligible magnitude compared to h 2 k 2 /2n when k itself 
is very small. 

The scattering length a is the limit of —A: -1 tan5fc for 
vanishing k [^|. The scattering length a(R c ) = a c is the 
corresponding limit of — k^ 1 tan Sk{R c ) and a c — > a as 
R c — ► oo. Calogero jg derived the Riccati equation in 
R c satisfied by a(R c ) or a c . We simplify the notation by 
dropping the subscript c and note that Calogero's def- 
inition of a(R) is opposite in sign to that used in the 
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effective range expansion 
equation is 



We see that the Riccati 



^ = [R-a(R)fv(R). 



(3) 



Hence we can find the scattering iength by solving Eq. 
over the interval [0, R c ] with initial condition a(0) = 
where R c is chosen sufficiently large. 



A. Change of Variables 



We make the change of variables |2f 



a{R) = tan 6>(i?), 



(4) 



and 



III. NUMERICAL METHODS 

The numerical solution of Eq. Q is non-trivial. The 
function a(R) contains poles that correspond to the 
bound states supported by the potential as we illustrate 
in figure 1 for the test case described in section 3.3. 
Finite difference methods such as the Runge-Kutta 
method are unsuitable. Two possible ways to solve 
Eq. © are to change its variables, and to solve a closely 
related equation by the log-derivative method. 
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FIG. 1: The accumulated scattering length, a(R), as a func- 
tion of interatomic separation R. 



R = tan <j)(R), 



(5) 



to obtain the equation 



sec 4 (j>sm 2 [0(4) - <f>)] V[tan(^)], (6) 

which we solve over a range [0, <f> c ] by the Runge-Kutta 
method 0|. The value of <j) c corresponds to R c and is 
close to 7r/2. 



B. Log-derivative Method 

In Eq. J3J) we substitute u(R) = [R — a(i?)] _1 to obtain 
the equation 



du(R) 
dR 



+ u 2 {R) - V(R) = 



(7) 



which is the Riccati equation for the log -derivative of the 
radial wavefunction at zero energy [lj, . The substi- 
tution is also inferred by the asymptotic form, R — a, of 
the radial wavefunction at large R. The lo g-d erivative 
method can be evaluated by the propagator E3 



hu{R+h)+ - h 2 V(R+h) 



1 + hu(R - h) + - h 2 V(R - h) 



h 2 V{R) 



h 2 V(R) 
(8)" 



with the initial condition that u(0) is very large, where 
h is the step length; the local truncation error is 0(h 6 ) 
[ft| . The success of the log-derivative method in treating 
the poles of u(R) is attributable to its derivation from fi- 
nite difference approximations to the wavefunction and 
its derivative [l6|, i it is an example of a symplec- 

tic integrator in which the distance coordinate and its 
canonical partner are advanced simultaneously 



C. Numerical examples 



We tested the numerical methods of sections 3.1 and 
3.2 by evaluating the scattering length for a pair of cae- 
sium atoms, each of mass 2.422 x 10 5 atomic units (a.u) 
interacting via the model potential described by Gribakin 
and Flambaum [l9| 
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V(R) = aR? exp(-jR)-(C 6 R~ 6 + C 8 



where 



C W R~ W ) f(R) 
(9) 



f(R) = H(R-R')+H(R'-R) cxp [-{R'/R - l) 2 ] (10) 

with a = 0.0008, (3 = 5.53, 7 = 1.072, C 6 = 7020, 
C* 8 = 1.1 x 10 6 , do = 1.7 x 10 s and R' = 23.165, all 
in atomic units. As is the case in Numerov's method 
the local truncation error in the log-derivative method 
is approximately proportional to h 6 \V(R)\ 3 |20|. To 
economise on computation and reduce possible trunca- 
tion error we doubled the step length whenever the local 
truncation error had reduced by approximately a factor 
of 2 0, 0| when using the log-derivative method; the 
propagator in Eq. (JSJ) is easily modified to accomodate 
the change in step length. When solving the first order 
equation, Eq. 10, we used the Runge-Kutta method 
with a self-adjusting step length. Both methods yielded 
the same scattering length as expected which was accu- 
rate to 7 significant figures with R c — 40000 bohr in 
agreement with the findings of Marinescu [ll( . The vari- 
able phase function is illustrated by the curve for a c in 
figure 2; the last pole is located at around 100 bohr and 
for R >100 bohr the function is smooth and approaches 
its limit ultimately from above. 



IV. CORRECTIONS TO THE SCATTERING 
LENGTH 

Eq. © illustrates clearly the effect of neglecting the 
potential for R> R c since the equation can be recast as 



by parts, replacing da/dR by Eq. © whenever it occurs, 
thus increasing the order each time. We find 

poo 

E c = EP + 2 [R- a(R)} 2 V(R) 

x {[R-a(R)]W(R) -X(R)}dR, (13) 
where eP is the first order correction 

E^ = -(R c - a c ) 2 W c + 2(R C - a c )X c - 2Y C (14) 

where the subscripts denote that the quantities are eval- 
uated at R = R c and 



W(R) = J V(R)dR, 



X(R) 



W{R)dR, 



(15) 



(16) 



and 



Y(R) = / X(R)dR. 



(17) 



When examining the corrections at R = R c we assume 
that the potential is attractive for R > R c . The cor- 
rection E c is clearly negative so that a c approaches a 
from above as illustrated in figure 2. If the potential is 
a negative series of inverse powers then V c < 0, W c > 0, 
X c < and Y c > 0; if also either a(R) < and R > R c 
or R > a(R) > R c then the integrand in the last part of 



/ [R - a(R)} 2 V(R)dR+ / [R- a{R)} 2 V(R)dR. Eq. is negative and E c < eP . 
Jo Jr c 



Hence 



The first part of this equation is the scattering length 
a c = a(i? c ) calculated by any numerical method in which 
the computation is stopped at separation R c . The second 
part is the contribution made by the long range interac- 
tion; the choice of R c — 40000 bohr in the calculations 
of section 3.3 was needed to ensure that the second part 
was sufficiently small to give seven figure accuracy in the 
scattering length. 



The First Order Correction 



E c < ac + EP, 



(18) 



and therefore the first order corrected scattering length 
approaches a from above. This is illustrated in figure 
2 by the curve a( u ). The first order correction EP is 
seen from Eq. p4|l to be negative. Hence the first order 
corrected scattering length, ci c + eP , is a better approx- 



imation than is a c 
bounds. 



Both a r and a r 



E 



(i) 



are upper 



The correction is 



poo 

E c = E(R C ) = [R- a(R)} 2 V(R)dR. 



B. Higher Order Corrections 
We rewrite Eq. (fT5|) as 



(12) 



We wish to extract the term that is of first order in the 
potential strength. We evaluate the quadrature, Eq. (|12|l . 



E c = EP + £ c 



(19) 



where 
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£ c = 2 / [R- a{R)} 2 V(R) {[R - a(R)} W(R) - X(R)}dR 



(20) 



r 



represents higher order corrections. It is desirable to ob- 
tain a lower bound to the scattering length. We do this 
by partially taking account of the higher order correc- 
tions to obtain an expression that is correct to at least 



first order but also includes sufficient higher order terms 
to be a lower bound. We integrate Eq. I|20|) by parts to 
find 



J 



poo 

£ c = £ { c 2) + [R- a(R)] V{R) {[R - a(R)} X(R) - 2Y(R)} dR, 



(2) 

where £ c takes account of some, but not all, of the second order terms and is given by 



(21) 



J 



£ {2) = -(R c - a c fW 2 + 3(i? c - a c ) 2 W c X c - 2(R C - a c ){X 2 c + W C Y C ) + 2X C Y C . 

I 



(22) 



The integrand in Eq. Ij21|l is positive and hence £ c > £c 2 \ 
Therefore 



where E^ +£^ is taken as the leading term (the ratio) 
on the right hand side of Eq. I|24|) . 



a>a c + EP + £^ 



(23) 



and the right hand side of inequality Eq. <2.'}|) is a lower 
bound. The expression E, 
Eqs. (O and (H), as 



(i) 



(2) 

£c can be rewritten from 



E, 



(i) 



-2Y C + 2{R C - a c )X c - (R c - a c ) 2 W c 
1 + X C -{R C - a c )W c 



where + ... denotes terms of order three and more in 
the potential strength. The terms Ec + £c^ constitute 
a guaranteed lower bound to the correction E c . The first 
term (the ratio) on the right hand side of Eq. Q24JI is 
not a guaranteed lower bound but, on the assumption 
that R c is sufficiently large that the higher order terms 
are negligible it is, in practice, a lower bound. We have 
achieved bounds on the scattering length; from Eq. (IT^) 
we have an upper bound 



o'°)=o c + ^ 1 ), (25) 
and from Eq. I|24|l we have a lower bound 



a^=a c + E^+£M 



(26) 
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FIG. 2: The accumulated scattering length, a c , and corrected 
scattering lengths, a' L ' and a< u >, as functions of interatomic 
separation R. 

Marinescu derived a series of corrections to the 
scattering length for differing orders of the potential 
strength. He derived corrections from an iterative so- 
lution to a sequence of differential equations, one correc- 



5 



tion for each order of the potential strength. On substi- 
tuting the functions defined by Eqs. JTSJ, and l(T7|) 
into Marinescu's first order correction we find that it is 
the same as the correction given by Eq. H26|) . Marinescu 
showed that his higher order corrections alternate in sign 
giving upper and lower bounds. Our equations l|25l) and 
(|26[) provide bounds very simply; while they are guaran- 
teed correct only to the first order they provide good, use- 
ful corrections for calculations with a moderately small 
value of R c . In the model caesium collision problem dis- 
cussed in section 3.3, range of 1250 bohr proved sufficient 
to yield seven figure accuracy. This range agrees with 
that found by Marinescu. The higher order corrections 
are complicated. The curve in figure 2 illustrates 
how the lower bound corrections leads to a scattering 
length that approaches a from below. 



Eq. (|28() agrees with our previous study [j| and, for n — 6, 
with the analysis of Hinckclmann and Spruch 10]. The 
correction to be made to a^ v \ given by Eq. (|2C)|> . is 
asymptotically, 



1 ■)) 



(n - 2)(2rc- 5) : 



and the correction for is 



#10 



~2n-4 



(rc - 2)(rc - 3)(2rc - 5) ~ 2(rc - 3) ' 



(31) 



(32) 



From Eq. 1\','>2\i we can construct a better approximation. 
In these asymptotic conditions it is 



C. Inverse Power Series Potentials 

For the calculations shown in figures 1 and 2, we used 
the asymptotic potential 



, _ a (u) + 2 (n - 3) a W 
Cl 2n-5 

We illustrate this in figure 3 for n=6. 



(33) 



V(R) = -C e R- 6 - C 8 iT 8 - C W R- 



(27) 



with the coefficients shown in section 3.3. With such 
inverse power potentials Eqs. i|14|) . (|24|l and l|2tj[l assume 
simple forms. 

When the term is R~ 6 is considered alone we find that 
Eq. i|14|) yields the expression that Szmytkowski ^3] de- 
rived from analysis of the asymptotic wavefunctions. In 
his analysis Szmytkowski gave also an expression for the 
scattering length in the form of a fraction involving Bessel 
functions. By replacing the Bessel functions by their 
asymptotic expansions for small arguments we find that 
Szmytkowski's ratio agrees with the lower bound in 
Eq. BB>. 

When a term —C n R~ n dominates the potential the 
upper and lower bounds are 



and 



5c - <5" 2 



1 



71 — 3 



2a c 
rc-2 



a 2 



n- 1 



(28) 
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FIG. 3: The corrected scattering lengths a' , aS > and a< u \ as 
functions of interatomic separation R. 



~n-2 



= a c - 



2a c 



n — 3 n — 2 n — 1 



1 - a„ 



rc — 2 rc — 1 



, (29) 



where the tilde notation indicates that lengths are re- 
placed by the dimensioneless quantities obtained by di- 
viding by i? c , and a n is a length characteristic of the 
potential given by 



V h 2 



l/(n-2) 



(30) 



V. SECULAR PERTURBATION THEORY 

Secular perturbation theory has been applied to mo- 
tion in the long range polarisation potential —C^R^^ by 
Cavagnero He showed how the scattering length is 
altered by long range interactions in accord with the pre- 
dictions of Spruch et al [2^] . We indicate below how sec- 
ular perturbation theory yields the lower bound a 1 -^ for 
the scattering length for a van der Waals —CqR~ 6 poten- 
tial. The result is easily generalised for a potential such 
as that in Eq. . 
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We are interested in the solution of Eq. (JTJ, with 
V(R) — —CqR~ 6 , for small values of the wavenumbcr 
k and large val ues of the separation R. Substituting 
R^yfcR) = y/(k/z)M(z), where z = kR, in Eq. Q 
we obtain 



z 2 - 



2[i k 4 C 6 



M(z) = - a 4 



M(z). 
(34) 

If the right hand side of Eq. (|34|1 were absent, M(z) would 
be a Bessel function of order 1/2. We replace M(z) by 
the series 



M(z) = ^ CnJv+2n(z), (35) 
n=— oo 

and, on the right hand side of Eq. $Q. we substitute |23l ] 



n+2 



l Jy+2n(z) = V m , n (v) Jv+2m(z), (36) 



m— n— 2 



where V m ,n(v) are expansion coefficients, to find linear 
algebraic equations for c„ 



n+2 

E 



9 1 

(1/ + 2m)" - - 



<Vn + ^ Kn,n(^) j C„ = 0, 

(37) 

where we have truncated the sum in expansion at 
n = ±n maa: . Eq. I|37|l leads to a determinantal equation 
satisfied by v. We include a term in k in 1/ and write 

M 



bk" 



(38) 



and, with this given value of u, we calculate the coef- 
ficients c n from Eq. Ij37(l after arbitrarily setting Co to 
unity. The equation with n = provides a new value for 



low energy phase shifts for the model caesium scattering 
problem were determined with a matching radius of only 
50 bohr ||. 

The secular perturbation expansion for the scatter- 
ing length appropriate for the potential in Eq. <|27|) is 
identical to that of Marinescu to the first two or- 
ders in the secular expansion. When the leading term 
dominates, the secular expansion yields the correction of 
Szmytkowski [l2l| . 

VI. CONCLUSION 



The variable phase method for potential scattering 
yields a simple first order differential equation satisfied by 
a function that tends in the limit of infinite interatomic 
separation to the scattering length 6] . The equation can 
be solved by specially adapted numerical methods, thus 
providing a method to compute the scattering length. 
This, and all other computations, being finite necessar- 
ily approximate the scattering length by its value com- 
puted at some maximum separation R c . The differential 
equation is readily employed to derive corrections aris- 
ing from long range interactions over the range [i? c ,oo]. 
We derived two such corrections accurate to at least first 
order which provide upper and lower computed bounds 
to the scattering length. In a model collision problem 
these corrections provide a 40-fold reduction in the value 
of R c that is necessary to let us obtain an accurate scat- 
tering length. The corrections are very simple and we 
suggest that use of the first order correction with a mod- 
erate value of R c is sufficient. Other methods involve 
a more complete but more complicated correction, such 
as the 50 th order secular perturbation theory 0, or 
higher order expansions of the solution of the zero-energy 
Schrodinger equation applied to a calculation with 
a much smaller value of R c . However, with techniques 
such as interval doubling 0, [2l| there is little advantage 
in having R c very small. We note that the first order 
differential equation for the effective range 6] , while not 
useful numerically, does yield the first order long range 
correction derived previously 0, llOfl . 
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